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Abstract 



Earlier study of quark-hadron phase transition in the Ginzberg-Landau the- 
ory is reexamined in the Ising model, so that spatial fluctuations during the 
transition can be taken into account. Although the dimension of the physical 
system is 2, as will be argued, both d = 2 and d = 4 Ising systems are stud- 
ied, the latter being theoretically closer to the Ginzberg-Landau theory. The 
normalized factorial moments F q are used to quantify multiplicity fluctuations, 
O ! an d the scaling exponent v is used to characterize the scaling properties. It 

is found by simulation on the Ising lattice that v becomes a function of the 
temperature T near T c . The average value of v over a range of T < T c agrees 
with the value of 1.3 derived analytically from the Ginzberg-Landau theory. 
Thus the implications of the mean-field theory are not invalidated by either the 
introduction of spatial fluctuations or the restriction to a 2D system. 



1 Introduction 
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It has recently been found in theoretical investigations that in quark-hadron phase 
transition (PT) there are hadronic observables that exhibit scaling behavior 0-0. It 
is based on the Ginzburg-Landau (GL) description of second-order PT. Subsequently, 
the consideration has been extended to first-order PT 0-0. Besides the unique fea- 
tures of the QCD dynamics, what distinguishes the problem from the standard criti- 
cal phenomena such as ferromagnetism is the lack of control of the temperature in a 
heavy-ion collision and the availability of only the hadrons as observables. Together 
with the complications in the nonperturbative regime of QCD, these difficulties pre- 
clude at this stage any observable predictions that follow from first principles. For 
that reason the Ginzburg-Landau theory has been used as the framework for quark- 
hadron PT, in which quantitative predictions on the hadronic observables can be 
derived and compared with experiment H [| . However, it is a mean-field theory that 
is accurate only in high dimensions. In the Ising model the upper critical dimension 
is four, above which the critical exponents determined in the theory is exact. The 
physical problem in heavy-ion collision, however, has dimension d — 2, as we shall 
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argue below. It is therefore important to consider the Ising model in d = 2 for quark- 
hadron PT, and compare its scaling behavior, if existent, with that derived in the 
GL theory. It is also useful to consider the d — 4 case so that the scaling result can 
be checked against the mean- field result. In that way one can gain some insight on 
whatever discrepancy there may exist between the GL and d = 2 Ising results. Those 
are the problems that we shall investigate in this paper. 

The use of Ising model as a solvable problem to study particle production has 
been done by many authors in the past several years |10]-[2(|. The focus has been 



on the issue of intermittency pli , i.e., power-law behavior of the normalized factorial 
moments as functions of the bin size. Differences in their conclusions appear to 
depend largely on how a particle is defined in terms of the lattice spins and how the 
limit of large N (number of lattice sites) is taken. In those papers that are specifically 
concerned with the d — 1 Ising model, no PT is contemplated. In the GL formalism 
it has been found analytically that there is no intermittency UJ; nevertheless, there is 
a different scaling behavior, called F-scaling, for which a numerical scaling exponent 
v has been determined. We want to check the value of v in the Ising model. 

One of the shortcomings of the GL approach is that the phase transition described 
by it is spatially too smooth. In a heavy-ion collision one expects hadrons to be formed 
one at a time or in clusters at the surface of a quark-gluon plasma (if created) with 
large spatial fluctuations. The Ising model is well suited to improve on that aspect of 
the problem. We shall consider cells as basic units consisting of several lattice sites, 
and vary the size of bins that contain certain multiples of the basic cells. In that way 
the model problem simulates the real situation in a much better approximation. The 
fluctuations that we study cannot be analyzed analytically, so computer simulation 
will be essential in carrying out the investigation. Our study in d = 4 can help to 
elucidate the question whether the difference between the GL and d = 2 Ising results 
is due to the specific dimension chosen or to the modified definition of hadrons in the 
Ising system. 



2 Hadronization on the Ising lattice 

The conventional picture of a heavy-ion collision that creates a quark-gluon plasma is 
an expanding cylinder of quarks and gluons with temperature T high in the interior, 
lower on the surface. If the plasma undergoes a second-order PT to hadrons, then 
there is no mixed phase, and the surface of the cylinder is where the PT takes place 
at T < T c , the critical temperature. After the hadrons are formed, they leave the 
plasma and are decoupled from the system under consideration. Thus focusing on 
the cylinder surface only, we are dealing with a two-dimensional system of quarks 
and gluons, which has patches that hadronize but are immediately replaced by the 
same quark system underneath after the hadrons are removed. The hadronization 
patches fluctuate in location and in size, as time evolves. What we shall analyze 
is the spatial fluctuation at any given time. If time is integrated over, which is 
equivalent to collecting all the emitted hadrons from the beginning to the end of the 
hadronization process, then the fluctuations we seek would be smoothened out by the 
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overlap of hadrons produced at different times. Thus in order to see the fluctuations 
in the experimental data (apart from the complications arising from the possibility 
of hadron gas, which is a totally different issue, not to be considered here), it is 
necessary to make px cuts. That is based on the general feature that the average 
Pt decreases with increasing time due to the decrease of transverse pressure, as the 
longitudinal and transverse expansion progresses. So selecting a narrow range of pt 
has the effect of selecting a narrow time range of hadronization, during which the 
spatial fluctuations can be frozen in the data p2| . 

The above discussion forms the basis of our consideration of a two-dimensional 
system for the second-order quark-hadron PT. We unfold a finite section of the cylin- 
der at midrapidity and map it onto a square lattice. To treat the PT problem on 
this lattice in the framework of the Ising model requires a careful consideration of 



what constitutes a hadron. Various definitions have been used in Refs. [|10|-[|20|; the 
most common one associates spin up to a particle at a site. That may be satisfactory 
for particle production without PT, but is unsatisfactory for our purpose because 
for T > T c there should be no hadrons produced except for those produced due to 
fluctuating effects. Our definition is the result of a long chain of reasoning, which 
goes as follows. 

Our starting point is the Ginzburg-Landau (GL) theory of PT considered in Ref. 
[0. One may question the validity of that starting point for hadron production. 
The support for that position is to be found in the production of photons at the 
threshold of lasing, which is known to be a second-order PT describable by the GL 
theory Indeed an experiment in quantum optics has accurately verified the 



prediction in |l[ of the scaling exponent v = 1.304 po| . It is reasonable then, in the 
absence of a workable theory to describe the fluctuations of hadronic observables, to 
replace photons by hadrons in the GL theory. In the laser problem the rate of photon 
emission is 

p(t) = \m 2 , a) 

where <ft(t) is the complex eigenvalue of the annihilation operator that specifies the 
coherent state \(f>). The probability of detecting n photons in time r with detector 
efficiency rj is then 



P n {r) = J d<f>d(f)*P° r]£ p{t)dt 



B{4>) (2) 



where P° [rjn{r)\ is the Poisson distribution corresponding to the multiplicities in a 
pure coherent state \<f>) . The integration in (fj) is weighted by the probability B(<f)) 
that the system is in the state \4>). That probability is given by a Boltzmann factor, 
whose exponent for a homogeneous single-mode laser has exactly the same form as 
the GL free energy for second-order PT [p3]|-[p5[ . 

We now adapt this formalism for hadron production with the identification || 

p{z) = \<P{z)\ 2 (3) 
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as the hadron density at the spatial coordinate z defined in a ID space. The Poisson 
distribution is explicitly 



p *=hXL dzp{z) ) n ^ [-J A dzp{z 

and the hadron distribution as a result of PT is |3J 

P n = Z- 1 jv<pP»[<P}e- F ^ , 

where Z = J T)(pe~ F ^\ and F[(f>] is the GL free energy 



(4) 



(5) 



dz 



a\ct){z)\ 2 + b\(j){z)\ A + c\d(j)/dz\ 



(6) 



Although the gradient term in (||) is explicitly taken into account in M , it is generally 
not considered in [@-@. The integrations over z in (|j) and @ are over the area of a 
bin of size S 2 . Thus the fluctuations examined are variations from bin to bin without 
further details within each bin. We want to improve on that by use of the Ising model. 

In the GL theory of ferromagnetism [2fJ the free energy is as given in (Eft with the 
order parameter <p(z) being identified as the local mean magnetization field, which in 
the Ising model is given by 



4- 1/2 E 



(7) 



where Sj is the spin ± 1 at site j. The sum in (^) is extended over an area A e 
(surrounding z) smaller than the resolution of the detector, which in our case is the 
bin size 5 2 , but large enough to contain many sites. Hereafter, we shall call that area 
a cell, having the size A e = e 2 , where e will be specified below. Thus in connecting 
quark-hadron PT with the Ising model, using the laser problem and the GL theory 
as intermediaries, we should identify the hadron density (|3|) with the absolute square 
of (0) so that the hadron multiplicity in the ith cell is 



n,; 



E 



S 3 



(8) 



where A is a scale factor that relates the quark density of the plasma at T c to the 
lattice site density in the Ising model. A relates a physical space to a mathematical 
space, and there is no way to determine a priori what that scale factor is. Any 
observable consequences of our calculation should, of course, be independent of A. A 
Heavyside function has been appended in (]8|) for reasons to be explained below. 
The Hamiltonian of the Ising model without external field is 



H 



(9) 
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where the sum is over the nearest neighbors. In the simulation of the spin configu- 
rations the net magnetization mi = J2ieL 2 s i °f the whole lattice of size L 2 may be 
either positive or negative. If tul < 0, we reverse all spins, replacing s, by — Sj at 
every site in the lattice, so that we always have w,l > 0, essentially by definition. We 
then define hadron multiplicity in a cell by (^), in which the ^-function puts rij = 0, 
when the net spin in cell i is not positive definite. Thus at T > T c the disordered state 
tends to have n.j near zero; even if ^ for certain cells due to fluctuations, there 
should be nearly as many cells in a bin that have no hadrons. At T < T c the ordered 
state has mostly nonempty cells. These are the characteristics that are physically 
reasonable, and should be incorporated in the description of the heavy-ion collisions 
where hadrons are produced from the cylinder surface in patches of fluctuating posi- 
tions and sizes. It is an aspect of the hadronization process that cannot be addressed 
in the GL theory. 



3 Scaling behavior in the Ising model in d = 2 

Monte Carlo simulation of the Ising model by use of the Metropolis algorithm suffers 



greatly from the phenomenon of critical slowing down, when T is near T c [26]], i.e. the 
autocorrelation time in generating new configurations becomes very long. Previous 
numerical work on the study of intermittency in the Ising model [|1(], |12| must have 
encountered that difficulty We avoid the problem entirely by adopting the Wolff 
algorithm [|27|j , with which there is no dependence on the number of steps required to 
remove autocorrelation. 

We work with a 2D lattice of size L = 288 and consider cells of size e = 4. We use 
10 4 sweeps to set up an initial configuration, and then generate 10 4 configurations to 
calculate event averages. The average hadron multiplicity (n) e in a cell in units of A 
is shown as a function of T in Fig. 1. The sharp drop occurs at T = 2.3 (in units of 
J/fcs), a value slightly higher than the analytical value of T c , 2.27, for infinite lattice 
[ffijfl . Our determination of the precise value of T c will be given below according to 
a well-defined procedure. In Fig. 1, (n) e /\ approaches 256 as T — > 0, as it should, 
since when all spins are aligned, the total-spin squared in a cell is (4 x 4) 2 . In the 
large T limit (n) e is not zero because not every cell has equal numbers of spin-up and 
-down. However, the nonvanishing of (n) e is not important for our result, for we shall 
examine scaling behaviors only in the neighborhood of T c . 

The observable whose scaling properties are to be studied is the normalized fac- 
torial moments [^1|] 

F q = ( — — ~ -*~) " ' (™ s ~ + lUjA . (10) 

\ \ n t)h lv 

where n$ is the hadron multiplicity in a bin of size S 2 : 

N S 

ns = J2rii , (11) 

i=l 
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rii being the cell multiplicity given in (ffl), the sum being over all N$ cells in a bin, where 
Ng = (5/e) 2 . In (|l(]) (• • •), means (horizontal) average over all bins in a configuration, 
and (• • -) v means (vertical) average over all configurations. The number of bins M in 
the lattice is M = {L/5) 2 = (L/e) 2 /Ns. Intermittency refers to the scaling in the bin 
size or in M 

F q {M) oc (12) 

over a range of M. More precisely, we shall refer to it as M-scaling to be distinguished 
from F-scaling, which is the power-law behavior 

F q cc f£\ (13) 

when M is varied. Clearly, the validity of (|TJ) guarantees (0), but flT3"| ) may be valid 
even if (|2|) is not. The latter situation is, in fact, the case in the GL theory 0. 
Phenomenologically, there is support for F-scaling, when the evidence for M-scaling 



is weak, from data that have no relevance to quark-hadron PT [29f|. 

To determine the F q moments we generate a spin configuration at a given T by 
the Wolff algorithm and calculate the cell multiplicity rij for a wide range of A for all 
cells in the lattice, which is partitioned into {L/5) 2 bins. Since for any given A the 
bin multiplicity ns need not be an integer, we take the largest integer less than n$ and 
use it in (|T0|) for each bin. We then perform the horizontal averages in (|T0|) over all 
bins in the lattice, followed by averaging over 10 configurations at the same T. The 
dependence of F q (M) on A is shown in Fig. 2 for a number of representative values 
of q and M. As expected, there is nearly no dependence on A when A is large in the 
interval (0, 1), since when ns is large the factorial product n${n$ — 1) • • • (n«s — q + 1) 
is approximately equal to n\ and the A factors are canceled in the ratio. But when A 
is very small, then the effect on F q is significant. 

Since the critical point is independent of how we map the quark density to the 
lattice density, we first set A = 1 and determine the critical temperature T c . In Fig. 
3 we show the M dependence of F q (M) for A = 1 at three values of T: 2.20, 2.315, 
and 2.40, always in units of J/ks hereafter. It is evident that M scaling is good for 
T = 2.315, but not for the other two temperatures. The quality of straightness in the 
log-log plots can be made more transparent by comparing just the q = 6 moments in 
the same figure, as shown in Fig. 4. Since at T c clusters of hadronic patches of all sizes 
are formed, we expect F q to exhibit M scaling at T c as a characteristics of criticality. 
Thus we have a criterion for the determination of T c : F q (M) should show strict M 
scaling at T c . In Fig. 5 we show the x 2 of straightline fits of all F q (2 < q < 6) for a 
range of T in the neighborhood of small \ 2 and find that \ 2 is lowest at T = 2.315, 
which is therefore the T c in this problem. 

Having determined T c , we now investigate the A dependence of M scaling. In Fig. 
6 we show at T c the M-scaling behaviors for six values of A: 1, 0.1, 0.05, 0.03, 0.02, and 
0.01. For each value of q the set of curves appear lower for lower values of A. Clearly, 
M scaling is still good up to M ~ 100 even for the lowest A, but for M > 100, F q 
bends over and even decreases at high M and small A. This behavior is reasonable, 
because at small A, the cell multiplicities n« are small, so the bin multiplicities ns are 
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also small at high M; thus the factorial product in ( [TUP vanishes for many cells at 
high q, and F q receives nonvanishing contribution only from the rare cells and events 
that have large fluctuations. This phenomenon is not an artifact of the Ising model, 
but is the realistic situation in the data on multiparticle production (see, for example, 
the Workshop proceedings of |30| , |3T|). For this reason our use of the Ising model in 
the way we define hadron multiplicities constitutes a concrete improvement over the 
previous work on the Ising model on the one hand, and over the GL approach to 
quark-hadron PT on the other. 

From Fig. 6 we can determine the slopes (p q from the regions where M scaling is 
good. Even for A = 0.01, that can be done between M = 10 and 100. The results 
for A = 1.0 and 0.01 are shown in Fig. 7. Evidently, the dependence of tp q on q is 
remarkably linear, which we parameterize as 

¥>, = C(<Z-1) ■ (14) 

The slopes determined from Fig. 7 are 0.135 and 0.121, respectively, which can be 
summarized as ( = 0.128 ± 0.007. This result should be compared to that of the 



analytical work of Satz on the Ising model JTT|, where £ = 1/8 for d = 2. The 
agreement is excellent despite the fact that no reference to hadrons is made in Jll] . 
This gives us confidence in our procedure of relating quark-hadron PT to the Ising 
model. 

The above result is for T = T c only. As mentioned earlier, there is no M scaling 
at T 7^ T c . We now show that F scaling as defined in ( |I3| ) is nevertheless good for a 
small range of T < T c . Furthermore, the slope (3 q satisfies the formula 

(3 q = (q-lY (15) 



found in JIJ. But first, at T c the validity of (|T2| ) and (|T4D implies the validity of flilf ) 
with 

Pq = = q~ 1 • (16) 

Thus at T = T c , we have v = 1 exactly to the extent that ([TJ]) is exact. Of course, 
( |i~4"|) is only a fitted result that is very good, but not exact. 

A direct evaluation of F scaling can be accomplished by examining F q vs F2 in 
log- log plots. In Fig. 8 we show those plots for three values of T at A = 1: (a) 
T = 2.20, (b) T = T C = 2.315, and (c) T = 2.40. It is clear that F scaling is good 
for the first two cases, but not for the third. Hadronization at T > T c is only due 
to fluctuations, and we shall not consider that region any more below. On the other 
hand, PT can take place at T in a small region below T c , so our F-scaling analysis 
takes into account a wider range of hadrons produced in a PT, not just at T c . 

The values of the slopes (3 q in Fig. 8(a) and (b) are shown by square symbols 
in Fig. 9 as functions of q — 1 in the log-log plot . They are well fitted by straight 
lines, thus verifying the validity of (|T5|). The exponents v are found to be 1.068 for 
T = 2.315, and 1.565 for T = 2.20. These are the results for A = 1. 

To see the A dependence, we show in Fig. 10 the F-scaling plots for four values 
of T at and below T c and for various values of A. The range of A in which F-scaling 
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is manifest shrinks as T is moved farther away from T c . For those values of T and A 
where straight-line fits can reasonably be made in those figures, the slope parameters 
I3 q are determined. They all satisfy (|i"5|) very well. The case for A = 0.05 is shown 
by diamond symbols in Fig. 9. Thus F scaling, when it exists, can be represented by 
a value of the scaling exponent z/, which depends on T and A. In Fig. 11 is shown 
z/(T) for various A. When no value of v is given, it means that F-scaling is not good 
enough to provide unambiguous extraction of f3 q . For A = 1 the range of T where v 
can be determined is largest. Fig. 11 shows that v increases from 1.07 at T c to 1.67 
at T = 2.1. But for A = 0.02 only at T c can v be determined to be 1.02. The result 
at T c , i.e., ^(T c ) = 1.02 — 1.07, differs slightly from the earlier result of v — 1 that 
follows from flTBD . They are close enough to be compatible, considering the fact that 
both procedures involve fitting of approximate scaling curves in different variables. 

The physical situation corresponds to low values of A, since the hadron multiplic- 
ities in a cell should not be too large. Taking only A = 0.02 and 0.05, we show the 
corresponding values of v in Fig. 12. The variation is between v = 1.02 and 1.52 over 
a narrow range of T below T c . It straddles very well the value v = 1.304 (dashed 
line) determined in the GL theory |l[. We recall that in the GL calculation there is 
no dependence of v on the GL parameters a{T) and b (so long as T < T c ) and on 
the dimension d. It is an average result that is situated in the middle of the range of 
v{T) that we have now determined with the spatial fluctuations taken into account. 
The two results therefore give mutual support to each other. 

We have varied the cell size by reducing e from 4 to 3. All the results are found 
to be virtually unchanged. Indeed, a sensitive dependence of e would render our 
procedure of determining the scaling behavior of particle production unreliable and 
physically irrelevant. 

Since T is not subject to experimental control in heavy-ion collisions, we have 
no way to check the T dependence of v{T) empirically. Detected hadrons, even in 
a small pt interval may be created in a range of T below T c , so the observation of 
v — 1.3 for all hadrons emitted would be a confirmation of both the GL and Ising 
descriptions of PT. Our study here not only generalizes the PT problem to allow for 
spatial fluctuations of hadronization patches, but also demonstrates that smaller v 
means larger fluctuation, since it is known that the fluctuation is greatest at T c . 

4 Scaling behavior in d = 4 

Before going into d = 4, let us first review what we have learned so far. We have 
been trying to relate three physical processes and models: (a) quark-hadron PT, (b) 
GL theory of PT, and (c) Ising model in d — 2. From (b) is derived in the scaling 
behavior of F q , which is the measure appropriate for (a); however, being a mean- 
field theory (b) is too smooth to reveal spatial fluctuations of (a). Since the physical 
dimension for the PT is d — 2, (c) is used to examine the scaling behavior where the 
discreteness of the hadron multiplicity in a cell is emphasized. The result affirms the 
F-scaling behavior found in (b), but with the scaling exponent v dependent on T, 
which is not a feature of (b). At T c the value of v agrees with what can be inferred 
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from the analytic treatment of 2D Ising model |IT[, but when averaged over a small 



range of T < T c , it agrees with the analytic result derived on the basis of (b). Thus 
the disagreement between the original, strict results of (b) and (c) (i.e., u — 1.3 and 
1.0, respectively) is removed by our present treatment of hadron formation in the 
context of (a). However, there is a level of understanding still unattained thus far. 
In carrying out our calculation we have used the definition of hadrons on the lattice 
as described in Sec. 2 at the same time when the 2D Ising model is used to simulate 
second-order PT. The question then arises as to which one of these two aspects of 
the problem is more responsible for the result of the previous section, i.e., our way 
of counting hadrons or the two-dimensionality of the space. To answer this question 
we perform a similar calculation using the hadron definition of Sec. 2, but in the 
d = 4 Ising model, which is known to have critical indices in agreement with the GL 
theory. In so doing we can learn whether the T dependence of v is a result of discrete 
hadron counting (as would be the case if the behavior persists in d = 4) or of the 
two- dimensionality (if the behavior disappears in d = 4). 

In d = 4 we must be content with smaller lattice size. We have used L = 24 in 
each dimension, and considered a cell as being of volume e 4 , with e = 2. Thus a cell 
has 16 sites, as in the case of d — 2 where e = 4. Our bins of lattice volume <5 4 can 
contain Ng = {S/e) 4 cells, which can vary from Ng = 1 to 12 4 cells. 

As before, we use (8) to define the hadron multiplicity on the 4D lattice, except, 
of course, that the sum is over j e V e (i), all sites of a 4D cell at i. Since a 4D lattice 
is far from the reality of a quark-gluon system, there is no precise meaning to the 
value of Hi calculated by (8), apart from the physical interpretation that it plays the 
role of multiplicity in 4D. The scale factor A should, nevertheless, be present as in 
the 2D case to make the present calculation as close to the previous one as possible, 
so A can be varied over the same range. As in 2D, we would not want to think that 
a cell can contain as many as 256 particles. 

We continue to use (9) as the Hamiltonian, but in 4D, and to employ the Wolff 
algorithm for simulation. Since we have been unable to find in the literature an 
analytical determination of T c for the 4D Ising model, we have to search for the 
critical point without any guidance except by the procedure used in the previous 
section. 

In Fig. 13 we show the average cell multiplicity in units of A as a function of T. 
It is the result of averaging over 3 x 10 4 configurations after initial runs of 3 x 10 4 
sweeps. Evidently, the transition temperature is not as distinct as in Fig.l for d = 2. 
To determine T c precisely we have attempted to find the temperature at which there 
is M scaling as in Fig. 3(b). Since T c is independent of A, we set A = 1 during the 
search. Unfortunately, we have not been able to find any T at which F q exhibits 
strict power-law behavior in M, which is now M = (L/5) A = (L/e) A /Ns. Perhaps it is 
because in 4D the value of L at 24 is not large enough to allow the full manifestation of 
long-range correlation. Fig. 14 is an example of the M dependence of F q . That figure 
is actually for T = T c , the value of which is to be determined later. At neighboring 
values of T the behaviors of F q vs M are similar, so we omit their presentation. 

For any set of F q at a given T we can plot F q vs F%. Fig. 15 shows such a plot for 
the F q in Fig. 14. Fitting all such plots by the F-scaling formula (13) we can determine 
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the scaling exponent v. Fig. 15 shows the dependence of v on T. Remarkably there 
is a dip at T = 6.7. Mindful of a similar dip in Fig. 11, we declare the critical point 
to be at T c = 6.7. In the absence of any other method that we can apply for a precise 
determination of T c , the dip in v{T) emerges as the only feasible method. Knowing 
T c allows us to focus on the behavior of the scaling exponent in its neighborhood. 

Returning to Fig. 14 we can examine the M dependences for various values of A at 
T c . Note that there is essentially no difference between A = 1.0 and 0.1, just as it is in 
Fig. 6 for d = 2. But for A = 0.01 all F q fall below 1 for M > 10 4 and are not plotted 
in Fig. 14. For the range of M shown for A = 0.01, which corresponds to 4 < 5 < 12, 
the dotted lines almost show M scaling. In the F-scaling plot of Fig. 15, there appears 
to be no dependence on A for 0.01 < A < 1.0 up to the second highest values for F q . 
The last and highest points for log-F 2 > 0.6 are for 5 = e = 2 (M = 2 x 10 4 ) and for 
A = 1.0 and 0.1 only. The large gap from the second highest ones (5 = 2e) cannot 
be reduced because the bin size must jump discretely from two cells to one in each 
dimension. The values of v are determined by straight-line fits of the portion of Fig. 
15 corresponding to 5 > 4, i.e., excluding the highest points. 

The result on z^(T), shown in Fig. 16, exhibits some A dependence at T c . For 
the case A = 0.01, which is perhaps closest to the physical problem with the lowest 
(n) per cell, we have v m 1.3, in agreement with the result in GL theory [|IJ. That 
is gratifying, since this is our only way of checking that the scaling behaviors of 4D 
Ising model and of the GL theory agree. The fact that v is distinctly different from 
1.0 implies that the scaling exponent is not universal. As in the 2D case, F scaling is 
valid only within a very small range of T around T c , when A = 0.01. 

When A > 0.1, we can investigate the T dependence of v(T). Fig. 16 shows that 
i/(T) increases as T decreases from T c , exhibiting essentially the same features as in 
Fig. 11. This result answers the question we have asked toward the end of the first 
paragraph of this section. The T dependence of v{T) is independent of the dimension 
d. Since the analytic result from the GL theory is that v — 1.3 independent of T ||T]|, 
we now see that the T dependence of v{T) is a consequence of the spatial fluctuations, 
which have been simulated on the lattice, but are lacking in the GL theory, and has 
little to do with the dimension of the system. 

5 Conclusion 

By working with the Ising model we have learned a number of things about second- 
order PT that are relevant to the quark-hadron problem. First of all, because the 
physical system of interest is two dimensional, it is the result of the 2D Ising model 
that we should concentrate on. Using the definition of hadrons on the lattice that 
we have constructed on the basis that hadron density is related to the square of the 
order parameter (or the total magnetization in a cell), we have been able to introduce 
spatial fluctuations during the hadronization process, a feature that is missing in 
the GL description. In so doing we are faced with a dilemma that the result of our 
simulation resolves. 

The dilemma is that the scaling behavior of 2D Ising model is known from analyt- 



10 



ical studies [IT] to imply v = 1.0, whereas the GL theory yields the analytical result 
of v = 1.3 [lj. The resolution is that when discrete hadron counting is achieved by 
lattice simulation, a temperature dependence is introduced into v{T). At T = T c , 
v{T c ) is indeed approximately 1.0, but at T < T c , v{T) is larger. The average is 
approximately 1.3, as shown in Fig. 12. In the GL theory v is independent of T 
so long as T < T c , so its result u — 1.3 may be regarded as an average over T in 
agreement with our simulated result. In the physical problem of quark-hadron PT 
in heavy-ion collisions, the hadronization temperature is not measurable. The ob- 
served hadron multiplicity fluctuations may well reflect hadronization occurring at a 
range of T < T c . Thus the average v = 1.3 is still the relevant scaling exponent in 
phenomenology. 

The T dependence of v{T) is not a special feature of d = 2. By simulating in d = 4 
we have found essentially the same T dependence. Moreover, at T c in 4D the value 
of ^(T c ) is approximately 1.3 rather than 1.0, reconfirming the known fact that the 
GL theory and the Ising model agree at d = 4 or higher. This gives us confidence in 
our choice of definition of hadron multiplicity in the Ising model and in the physical 
implications of the results obtained by simulation in 2D. 

The basic understanding we have learned is that when the shortcoming of the GL 
theory is amended by taking into account the spatial fluctuations of hadronization, 
the predicted result of v = 1.3 is unchanged, except for the realization that it is a 
mean value after averaging over a range of T below T c . 

It is also important to recognize at a more general level that in our attempt to 
relate quark-hadron physics to statistical physics we have found the effectiveness of 
using the factorial moments F q , the significance of F-scaling, and the amazing validity 
of (3 q = (q— l) v in all circumstances, with the consequence that the scaling exponent v 
has emerged as an effective parameter characteristic of PT, in lieu of the usual critical 
indices that rely on T being under experimental control. It is therefore suggestive that 
our method of describing the properties of PT may find applications in other areas of 
science where discrete multiplicities are the only observables at one's disposal. Such 
an application has recently been found in evolutionary biology 
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Figure Captions 

Fig. 1 Average multiplicity in a cell ( in units of A) as a function of temperature 
for d = 2 lattice. 

Fig. 2 Normalized factorial moments as functions of the scales factor A for q = 3 
and 6 and for three values of M. 

Fig. 3 F q vs M at three temperatures. 

Fig. 4 Same as in Fig. 3 but for q = 6 only. 

Fig. 5 x 2 °f straight-line fits of logF q vs logM for T in the neighborhood of 2.315. 

Fig. 6 M-scaling plots of F q at T c for various values of A. In each cluster for a 
fixed q lower curves have lower A. 

Fig. 7 Intermitting indices (f q at T c for A = 1.0 and 0.01. 

Fig. 8 F-scaling plots at three temperatures. 

Fig. 9 The slope f3 q in log- log plot for different T and A. 

Fig. 10 F-scaling plots for different T and A. 

Fig. 11 Scaling exponent v(T) for different A. 

Fig. 12 z/(T) at small A compared to the analytical value v = 1.304 from [1]. 

Fig. 13 Average multiplicity in a cell ( in units of A) as a function of T for d = 4 
lattice. 

Fig. 14 M-scaling plots for different A at T c in 4D. 
Fig. 15 F-scaling plots for different A at T c in 4D. 
Fig. 16 v(T) near T c for different A in 4D. 
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